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ABSTRACT 

It has long been speculated that the nature of the hard X-ray corona may be an important second driver of 
black hole state transitions, in addition to the mass accretion rate through the disk. However, a clear physical 
picture of coronal changes has not yet emerged. We present results from a systematic analysis of Rossi X- 
ray Timing Explorer observations of the stellar mass black hole binary XTE J 1650-500. All spectra with 
significant hard X-ray detections were fit using a self-consistent, relativistically-blurred disk reflection model 
suited to high ionization regimes. Importantly, we find evidence that both the spectral and timing properties of 
black hole states may be partially driven by the height of the X-ray corona above the disk, and related changes in 
how gravitational light bending affects the corona-disk interaction. Specifically, the evolution of the power-law, 
thermal disk, and relativistically-convolved reflection components in our spectral analysis indicate that: (1) the 
disk inner radius remains constant at r- ln ~ 1.65 ± 0.08 GM/c 2 (consistent with values found for the ISCO 
of XTE J 1650-500 in other works) throughout the transition from the brighter phases of the low-hard state to 
the intermediate states (both the hard-intermediate and soft-intermediate), through to the soft state and back; 
(2) the ratio between the observed reflected X-ray flux and power-law continuum (the "reflection fraction", R) 
increases sharply at the transition between the hard-intermediate and soft-intermediate states ("cannonball" jets 
are sometimes launched at this transition); (3) both the frequency and coherence of the high-frequency quasi- 
periodic oscillations (QPOs) observed in XTE J 1650-500 increase with R. We discuss our results in terms of 
black hole states and the nature of black hole accretion flows across the mass scale. 

Subject headings: X-rays: accretion - accretion disks - black hole physics - line: profiles - relativistic pro- 
cesses - X-ray binaries - X-ray Individuals: XTE J1650-500 



1. INTRODUCTION 

The majority of stellar-mass black holes residing in low 
mass binary systems spend most of their lives in quiescence. 
This explains the empirical fact that out of tens of thousands 
of such systems predicted to exist throughout our Galaxy 
(e.g. Yungelson et al. 2006), only about 50 have been dis- 
covered. For a number of such systems, it is well established 
that the outburst - which resulted in their discovery - evolves 
through a number of spectral states characterised by the rela- 
tive strength of their thermal and non-thermal X-ray emission 
and with possible differences in the accretion geometry and 
reflection attributes. These active states can be roughly sepa- 
rated into four semi-distinct states which are phenomenologi- 
cally described below and extensively discussed in Remillard 
& McClintock (2006) and Belloni (2010). At the outset of the 
outburst, the system goes through what has been traditionally 
dubbed the low-hard state (LHS) where the X-ray spectrum 
is dominated by a non-thermal component often simply de- 
scribed by a power-law (photon index T between ~ 1.4 — 2) 
spectrum with relatively low luminosity (~ 0.05iEdd and and 
exponential cut-off at ~ 100 keV). The energy spectrum in 
the LHS peaks near ~ 100 keV and often we also see the weak 
presence of a thermal component (contributing < 20% of the 
total 2-20keVflux) with a temperature below <~ 0.5 keV pro- 
duced by the accretion disk (see e.g. Reynolds & Miller 201 1; 
Reis et al. 2010, and references therein). 
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As the luminosity increases, the spectrum moves through 
the intermediate state (IS) where the 2-10 keV flux is typi- 
cally a factor of ~ 4 times higher than that of the LHS. Here, 
the soft (T = 2 — 3) power-law tail coexists with a strong ther- 
mal component. Recently, the intermediate state has begun to 
be subdivided into an early Hard-Intermediate State (HIS) and 
a later Soft-Intermediate state (SIS) just prior to a transition 
into the canonical High-Soft or thermal state, where the X- 
ray flux is dominated (> 75% of the total 2-20keVflux) by 
the thermal radiation from the inner accretion disk having an 
effective temperature of ~ 1 keV. In this final state of the out- 
burst, the system usually emits with luminosities > 0.1L Edd 
and the power-law component is both weak (less than 25% of 
the total 2-20 keV flux) and steep (T = 2 - 3). Following the 
HSS, the system often returns to the LHS and subsequently 
goes back to quiescence were it can remain indefinitely or in 
some cases for a few years before this cycle restarts. 

The hard X-ray emission predominant in the LHS has long 
been linked to inverse Compton scattering of the soft thermal 
disk photons by a population of hot (~ 10 9 K) electrons in 
a cloud of optically thin, ionised gas or "corona" surround- 
ing the inner parts of the accretion disk (Shapiro et al. 1976; 
Sunyaev & Titarchuk 1980). Under the common assump- 
tion that the radio emission observed to originate from stel- 
lar mass black holes, is directly related to the presence of a 
jet, it is believed that all such systems, either in the LHS or 
in transition, launch a collimated outflow (e.g. Fender 2001; 
Fender et al. 2004, 2009). The fact that these persistent jets 
are observed only in the LHS suggests that the jet is linked 
to the corona, with claims that the corona in the LHS is in- 
deed the launching point of persistent jets (see e.g. Markoff 
et al. 2005). The connection between the radio (jet) and X- 
ray flux for both stellar-mass and supermassive black holes 
(Merloni et al. 2003; Gallo et al. 2003; Falcke et al. 2004), 
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often referred to as "the fundamental plane" of black hole ac- 
cretion, suggests an intimate connection between the corona 
and radio-jets (see e.g. Miller et al. 2012). Whether state tran- 
sitions are driven by intrinsic changes in m, physical changes 
in the disk, disk-corona, radio jet or a combination of all these 
factors is a matter of much debate. 

1.1. Reprocessed X-rays: Reflection 

The existence of a hard X-ray source - the corona - also 
adds further complexities to the various spectral states. The 
reprocessing of these hard X-rays by the relatively cold ac- 
cretion disk in all active states results in a number of "reflec- 
tion features" consisting of discrete atomic features together 
with a "Compton-hump" peaking at approximately 30keV. 
The high fluorescent yield - and relatively high cosmic abun- 
dance - of iron often results in a particularly strong feature at 
^6—7 kcV (see e.g. Fabian & Ross 2010, for a recent review 
of "reflection" in black holes). 

The strong irradiation of the black hole accretion disk by 
the coronal photons likely causes the surface layers to be pho- 
toionised. Ross & Fabian (1993) investigated the effect of 
allowing the gas constituting the top layers of the accretion 
disk to ionise, and the authors went on to compute reflec- 
tion spectra for different ionization levels. A number of sim- 
ilar studies of reflection from ionised matter have been con- 
ducted since (Matt et al. 1993, 1996; Ross et al. 1999; Bal- 
lantyne et al. 2001; Garcia & Kallman 2010; Garcia et al. 
201 1). These studies demonstrate that the reflection spectrum 
expected from a black hole depends strongly on the level of 
ionization of the surface layers of the disk. This can be quan- 
tified for a constant density gas by the ionization parameter 



where L x is the ionising luminosity of the source, d is the 
distance between the disk and the source, and n is the density 
of the disk. Thus an increase in £, either by increasing the 
illuminating flux, decreasing the density or distance between 
the X-ray source and the disk, will cause the gas in the disk 
to become more ionised. Matt et al. (1993, 1996) split the 
behaviour of the reflection spectrum into four main regimes 
depending on the value of £. 

For low ionization parameter (£ < 100 erg cm s _1 ), the 
material is only weakly ionised and the reflection spectrum re- 
sembles that arising from "cold" matter, with a prominent iron 
line at 6.4 kcV, and strong absorption below ss 10 kcV. There 
is only a weak contribution from the backscattered continuum 
at « 6 keV and a weak iron K absorption edge at 7.1 keV. As 
the disk becomes more ionised (100 < £ < 500 erg cm s _1 ) 
the system reaches intermediate ionization range where Fe has 
lost all of its M-shell (n — 3) electrons and thus exists in the 
form of FeXVII-FeXXIII with a vacancy in the L-shell of the 
ion. Due to this vacancy, the L-shell can absorb the Ka line 
photons and thus effectively trap the escaping photon. This 
resonance trapping is only terminated when an Auger elec- 
tron is emitted. This second ionization regime is therefore 
characterised by a very weak iron line and a moderate iron 
absorption edge. 

As the gas becomes more ionised (500 < £ < 
5000 erg cm s _1 ) all low-Z metals are found in their hydro- 
genic form and the soft reflection spectrum has only weak 
spectral features. Iron is found mostly in its hydrogen or he- 
lium like forms (FeXXVI or FeXXV respectively) and due to 



the lack of at least 2 electrons in the L-shell (i.e. not having a 
full 2s sub-shell), Auger de-excitation cannot occur. The re- 
sult is strong emission from "hot" Ka FeXXV and FeXXVI 
at 6.67 and 6.97 keV respectively and the corresponding ab- 
sorption edges at approximately 8.85 and 9.28 kcV respec- 
tively. Finally, when £ 3> 5000 erg cm s _1 , the disk is highly 
ionised and there is a distinct absence of any atomic features. 

A further complication arises in the reflection spectra of 
stellar mass black holes due to the fact that in these systems 
the gas in the accretion disk is inherently X-ray "hot" meaning 
that \ow-Z metals can be fully ionised in the gas even before 
receiving any irradiation by the X-ray corona. To account for 
this extra "thermal ionization", Ross & Fabian (2007) per- 
formed self-consistent calculations of the reflection resulting 
from the illumination of the accretion disk by both a hard, 
powerlaw corona and thermal disk blackbody radiation. The 
authors compared the results of having the disk both in hydro- 
static equilibrium and under the assumption of a constant den- 
sity atmosphere, and found reasonably good agreement be- 
tween the two emergent spectra. Furthermore, Ross & Fabian 
(2007) also confirmed in stellar mass black holes the result 
that had been previously found for AGN in that the spectrum 
from a constant density disk is slightly diluted ("spectrally 
hardened") in comparison to that of a disk in hydrostatic equi- 
librium. The two reflection grids resulting from the work of 
Ross & Fabian (2005, 2007), will be used frequently through- 
out this work. 

1.2. General Relativistic Effects: Light bending 

Naively, assuming isotropic coronal emission, one would 
expect variations in the reflection component to be directly 
correlated to variations in the observed power-law continuum. 
That is, as the observed flux of the X-ray corona increases, 
so should the amount of reprocessed emission. However, in 
a number of instances it has been found that this is not the 
case, with the reflection component at times behaving in an 
anticorrelated manner (e.g. Rossi et al. 2005) or not varying at 
all despite large variations in the X-ray powerlaw continuum 
(e.g. Fabian et al. 2002; Fabian & Vaughan 2003; Miniutti 
et al. 2003; Ballantyne et al. 2003; Larsson et al. 2007). 

By virtue of its proximity to the black hole, the emission 
from the corona is naturally affected by general relativistic 
(GR) effects. Some of the radiation from the corona which 
would otherwise escape is gravitationally focused - "bent" 
- towards the accretion disk giving rise to enhanced reflec- 
tion and selectively decreasing the X-ray continuum at infin- 
ity. A number of studies (for instance Martocchia & Matt 
1996; Martocchia et al. 2000; Miniutti & Fabian 2004; Mini- 
utti et al. 2004) have investigated the effect of GR on a com- 
pact, centrally concentrated X-ray corona close to a black 
hole 5 . The "light-bending" model put forward by Miniutti 
& Fabian (2004) predicts a number of semi-distinct regimes 
affecting the variability of the reflection component compared 
to the X-ray continuum: 

Regime 1: When the corona is very close to the black hole (a 
few gravitational radii r g =GM/c 2 ), a large fraction of 
the radiation is bent onto the the accretion disk thus sig- 
nificantly reducing the amount of observed X-ray con- 

5 Observational evidence for such compact X-ray corona has recently 
come from microlensing results where the size of the X-emitting region has 
been shown to be of the order of ~ 10r g (Chartas et al. 2009; Dai et al. 2010; 
Chartas et al. 2012; Morgan et al. 2012). 
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tinuum and enhancing the reflection. A very steep emis- 
sivity profile is expected as the source is highly concen- 
trated in the inner regions and the reflection is expected 
to be a steep function of the continuum in a quasi-linear 
manner. 

Regime 2: When the central corona is slightly further from 
the black hole (at heights of ~ 10r g ), light bending 
causes the reflection component to vary significantly 
less than the X-ray continuum. The amount of light 
bent towards the black hole decreases as the corona 
moves further from the black hole and the X-ray con- 
tinuum increases. 

Finally, at heights 3> 20r g , light-bending becomes less im- 
portant and the observed continuum increases. In this man- 
ner, the presence of gravitational light-bending has been in- 
voked to explain the fact that Seyferts (e.g. Fabian et al. 2009, 
2012b) and XRBs (e.g. Reis et al. 2012) at times appears to 
be "reflection-dominated". Sources where the observed X-ray 
spectrum display a distinct lack (or comparably small amount) 
of hard, powerlaw-like continuum despite a large contribution 
of reflection are fully consistent with the first Regime detailed 
above. 

Of course, the model presented above is idealised in that all 
characteristics of the observed variabilities are assumed to be 
a result of variation in the height of an isotropic and compact 
corona with a fixed luminosity. Although intrinsic variation 
in the luminosity of the corona may well be present, it is un- 
likely that they could solely explain the behavior of the re- 
flected emission described above. Indeed, the clear presence 
of broad and skewed iron emission lines in a growing number 
of sources ranging from stellar mass black holes (Miller 2007; 
Miller et al. 2009b; Reis et al. 2009a; Steiner et al. 2011; 
Hiemstra et al. 2011) to AGNs (Tanaka et al. 1995; Nandra 
et al. 1997, 2007; Fabian et al. 2009; Brenneman et al. 2011) 
(and also to lesser extent neutron stars (Cackett et al. 2008, 
2010; di Salvo et al. 2005, 2009; Reis et al. 2009b)) strongly 
attest that general relativistic effects play an important role in 
producing the line profile, further supporting the notion that 
the corona and the inner disk are both in the inner regions 
around a black hole. 

1.3. XTEJ1650-500 

Amongst one of the first systems to provide observational 
evidence for the aforementioned effect of gravitational light 
bending around a black hole and indeed the first around a 
stellar mass black hole was XTE J 1650-500, which was dis- 
covered by RXTE on 2001 September 5 as it went into out- 
burst (Remillard 2001). Based on the spectrum obtained early 
in decay of the outburst by XMM-Newton, and more impor- 
tantly on the presence of a clearly broad iron emission line, 
Miller et al. (2002) were able to not only infer that the ob- 
ject in XTE J 1650-500 was indeed a black hole, but also that 
it was close to maximally rotating with a dimensionless spin 
parameter a s=y 0.998. 

By decomposing the hard X-ray continuum from the 
reflection component in three BeppoSAX observations of 
XTE J 1650-500, Miniutti et al. (2004) were able to show that 
the latter remained nearly constant despite a large change in 
the direct continuum, in a manner consistent with the predic- 
tions of light bending around a black hole. Optical obser- 
vations obtained after the system had returned to near quies- 
cence (Orosz et al. 2004) revealed a mass function f(M) — 



2.73 ± 0.56 M Q with a most likely mass of ~ 4M , and 
with it secured XTE J 1650-500 as genuine black hole binary 
system. Corbel et al. (2004) reported on the radio and X-ray 
observations of XTE J 1650-500 during the outburst. The au- 
thors find a clear drop by nearly an order of magnitude in the 
radio flux at the transition from the hard intermediate state 
(referred to as the intermediate state in that work) to the soft 
intermediate state (referred to as the steep power-law state in 
that work), and surprisingly they find residual radio emission 
during the often radio quiet disk-dominated soft state which 
they attributed to possible emission from previously ejected 
material interacting with the interstellar medium, rather than 
originating in the central source. 

A follow up study of RXTE data by Rossi et al. (2005) used 
the iron line as a proxy to the total reflection component and 
confirmed the plausibility of the light-bending scenario for the 
evolution of XTE J 1 650-500. Again using data obtained from 
the RXTE, Homan et al. (2003) reported on the discovery of a 
~ 250Hz QPO together with a number of less coherent vari- 
ability peaks at lower frequencies. By studying the spectral 
and timing evolution during the first <~ 80 days of the out- 
burst, the authors were able to define six periods (I- VI; in this 
work referred to as P1-P6) having somewhat distinct spectral 
and timing characteristics (see their Fig 1 and Table 1). A re- 
cent study involving XTE J 1650-500 has discussed the simi- 
larities between X-ray binaries and AGNs (Walton et al. 2012) 
and argued that both XTE J1650-500 and the active galaxy 
MCG-6-30-15 (Tanaka et al. 1995) must contain a rapidly 
rotating black hole, with the spin of XTE J 1650-500 having 
been formally constrained to 0.84 < a < 0.98. A further 
body of work based on RXTE observations and a variety of 
empirical models for the hard X-ray continuum (Yan & Wang 
2012) has concluded that the emission region, here referred to 
as the corona, decreases by a factor of <~ 23 in size during the 
transition from the hard to the soft state. 

1.4. This Work 

Since its discovery, XTE J 1650-500 has become one of the 
best studied black hole systems. However, the energy spectra 
of this system have either been studied in a great degree of 
detail using high-quality, single snapshot observations with 
XMM-Newton (i.e. Miller et al. 2002; Walton et al. 2012) or 
BeppoSAX (Miniutti et al. 2004) or using mostly phenomeno- 
logical and simple models in the study of long term evolu- 
tions with RXTE (i.e. Rossi et al. 2005; Dunn et al. 2010, 
201 1 ; Yan & Wang 2012). In this paper, we use the full RXTE 
archival data of the outburst to investigate the evolution of 
the direct power-law continuum, reflection and thermal disk 
components using, for the first time, a fully self-consistent 
prescription for the reflection component. 

In this manner, we combine the virtues of detailed analyses 
of single observations with the immense diagnostic power of 
multiple RXTE pointings. By being able to decouple the total 
reflection component (Fe-Ka emission line together with all 
other reflection signatures) from the illuminating continuum, 
we find that the transition from the hard-intermediate state to 
the soft-intermediate state is accompanied by a sharp increase 
in the strength of the reflected emission in comparison to the 
direct continuum. We interpret this increase in the reflection 
fraction as a sudden collapse of the corona as the system ap- 
proaches the thermal state, although we note that this may not 
be a unique interpretation. 

This paper is structured as follows: §2 briefly introduces 
the observations and details our various selection criteria. §3 
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FIG. 1 . — The 181 observations used were taken during the 2001/2002 outburst which is clearly seen in the All Sky Monitor (left). Since its discovery outburst, 
XTE J1650-500 has remained in quiescence. (Right:) PCA-PCU2 (top) and HEXTE-A (bottom) count rate during the time encompassing the outburst. The 
different colors marks the distinct spectral states as defined by Homan et al. (2003) based on the timing characteristics of the source. 



describes the base model and assumptions used throughout 
this work. The various results are presented in §4 and in §5 
we present a qualitative picture of a possible physical scenario 
that combines all our findings. 

2. OBSERVATION AND DATA REDUCTION 

We downloaded and analysed all 181 individual RXTE 
pointed observations of XTE J 1650-500 following its discov- 
ery. This gave a total of 307.4 ks PCA exposure which were 
reprocessed with the latest HEASOFTv6.12 tools. Over the 
last decade or so, RXTE has proven to be one of the most reli- 
ably calibrated X-ray observatories with one of the best doc- 
umented data reduction procedures. As such, we followed all 
the well established, standard reduction guides 6 . Given that it 
is the only PCU that is always on, as well as being the best 
calibrated of all units, we used only data from PCU-2. We 
chose to use the standard offset value of < 0.02 as well as 
elevation above the Earth's limb > 10°. Background spectra 
were produced using PCABACKEST using filter files created 
using XTEFILT. The latest PCA-history file was also used and 
throughout this work we use PCA data between 3 keV (ig- 
nored channels < 6) and 25 keV without additional system- 
atic errors as we are mostly interested in relative changes and 
the impact on the errors of the various parameters as well as 
in the \ 2 distribution shown in Fig. 6 is minimal with the only 
change being a systematic shift to lower values. 

For our analyses, we require both PCA and High Energy X- 
ray Timing Experiment (HEXTE) data (but see below). The 
HEXTE data were also reduced in the standard manner fol- 
lowing the procedures outlined in the RXTE guide. Back- 
ground files were generated using HXTBACK and the spectra 
were corrected for deadtime using HXTDEAD. The appropri- 
ate response and ancillary files were downloaded 7 and added 
as a header key using the FTOOL FPARKEY. HEXTE data 
were fit between 25-150 keV. 

Figure 1 (left) shows the 1-day averaged long term light 
curve as seen from the RXTE All Sky Monitor (ASM) where 
the 2001 outburst is clearly visible. The panels on the right 
shows the PCA-PCU2 count rate (top) and HEXTE-A (bot- 
tom) during the time roughly encompassing the outburst. The 
colours highlight the various states during the outburst and 
are directly mapped into the hardness intensity diagram (HID) 
shown in Fig. 2. In short, the first few observations caught 

6 Found at http://heasarc.nasa.gov/docs/xte/recipes/cook_book.html 

7 From ftp://legacy.gsfc.nasa.gov/xte/calib_data/hexte_files/DEFAULT/ 
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FIG. 2. — Hardness intensity diagram made using the 116 observations 
that remain after imposing a cut where both HEXTE Cluster A and B have 
source-background counts greater than 0. The vertical dashed lines marks 
approximate transitions between LHS - HIS- SIS -HSS. The large diamond 
symbols mark the position in the HID diagram for the representative spec- 
tra shown in Fig. 4. The color code is identical to Fig. 1 (right) and will 
remain the same throughout this work. We have scaled the size of the sym- 
bols for each observation to the value of the reflection fraction R, as deter- 
mined in § 4.1 and shown in Fig. 3. The size-legends from top-to-bottom are: 
R = 0.5, 0.75, 1,2,3, 4, 5. The falling branch of the LHS and the HSS have 
sizes corresponding to R = 0.5 as R is poorly constrained in these states. 

the source in a rising LHS which evolved to the HIS approxi- 
mately 3 days later, where it remained for ~ 15 days until the 
clear change to the SIS. We have further divided the HIS into 
early Period 1 and late Period 2, similar to the division made 
by Homan et al. (2003) based on the timing characteristics 
displayed. The source then remains in the SIS (P3) for ap- 
proximately 12 days before it makes the typical excursion to 
the disk dominated HSS (P4) where we see a clear drop in the 
HEXTE count rate in Fig. 1 . After ~ 40 days, the hard flux 
in XTE J 1650-500 sharply increases (P5 and P6) as the sys- 
tem returns to LHS and eventually goes back into quiescence 
where it remains up to the present time. 

The work presented in this paper is fully interpreted within 
the framework of reflection models (see e.g. Fabian & Ross 
2010, and §1.1). As such we are mainly interested in being 
able to determine the contribution of the reflected emission to 
the overall spectra. As mentioned in the introduction, reflec- 
tion is not limited to the iron line profile and in order to obtain 
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the best handle on the full reflection component we restricted 
the analysis that follows by requiring that both the HEXTE- 
A and the HEXTE-B units have background-subtracted count 
rates that are greater than zero. This results in 1 16 PC A point- 
ing observations totalling 217.9 ks. In doing so, we are effec- 
tively reducing our sensitivity to observations in very steep 
or disk dominated states, as well as to faint LHS, and most 
of the results presented here concern the intermediate states, 
which are more luminous at high energies (see Fig. 1; right). 
In order to investigate the effect of the HEXTE data on the 
results presented in this paper, we have also repeated all the 
fits described below using only the PCA data and confirmed 
that, although the results obtained do not strictly depend on 
the inclusion of the HEXTE data, owing to its relatively small 
statistical weight, the high energy data still provide a useful 
additional test of the continuum model (see e.g. Fig 4). 

In the following section we detail the spectral fits to the 
1 16 observations. All our work makes use of the X-ray fitting 
package XSPECvl2.7.0 (Arnaud 1996). Where uncertainties 
on models parameters are quoted, this refers to the 90 per cent 
confidence limit for the parameter of interest. 

3. THE MODEL 

Previous attempts to characterise the evolution of 
XTE J 1650-500 as observed with RXTE have relied on 
a purely phenomenological interpretation of the reflection 
continuum and features 8 and often employed a single GAUS- 
SIAN emission line with centroid energy fixed at 6.4 keV 
as expected from neutral Fe-Ka (i.e. Dunn et al. 2011), a 
combination of a similar GAUSSIAN together with a broad 
smeared edge (SMEDGE, Ebisawa et al. 1994) component 
(i.e. Yan & Wang 2012), or in the more physically appropriate 
application, a relativistic emission line (LAOR, Laor 1991) 
together with SMEDGE (Rossi et al. 2005). However, even 
in the latter example, the model was unphysical as the 
combination of SMEDGE with an emission line does not keep 
consistency between the depth of the edge - the number of 
absorbed photons - and the strength of the line. 

In the study of Dunn et al. (201 1), where the authors were 
primarily interested in the behavior of the thermal accretion 
disk, such a simplification is easily justified as small devia- 
tions at high energies are unlikely to have significant effects 
at low energies. Rossi et al. alluded at the importance of using 
reflection models but unfortunately were swayed away due to 
the highly time-consuming and computer-intensive task that 
this would present nearly 10 years ago. 

The importance of reflection from accretion disks can be 
directly measured by the sheer number of theoretical works 
that have been devoted to fully characterising its behavior 
(i.e. Lightman & White 1988; George & Fabian 1991; Matt 
et al. 1991; Ross & Fabian 1993; Zycki et al. 1994; Nayak- 
shin et al. 2000; Ballantyne et al. 2001; Ross & Fabian 2005, 
2007; Garcia & Kallman 2010; Garcia et al. 201 1). A widely 
used reflection code is that of Ross & Fabian (2005, RE- 
FLIONX), which provides a self-consistent treatment of the 

8 Contrast this phenomenological approach to the reflection component 
with the systematic testing of Comptonisation models for the power-law 
(corona) component by Nowak et al. (2002) on RXTE data of GX 339-4. 
Ideally one would strive to combine physical models for both the power-law 
and reflection (and of course the accretion disk which is known to not be a 
simple multicolour disk due to various relativistic effects) however this quest 
for fully relativistic disk and reflection together with physical prescription 
for the Comptonisation continuum is proving highly complicated even for 
single, dedicated efforts at snapshot observations and is beyond the scope of 
this paper. 



dominant atomic processes around a black hole and, given 
an input power-law continuum, outputs a reflection spectrum 
where both the "Compton-hump", emission and absorption 
features are all physically linked. 

The reflection spectrum is calculated in the local frame, 
therefore we employ the convolution model KDBLURf (Laor 
(1991); vastly optimised by Jeremy S. Sanders and employed 
in Fabian et al. 2012a and Reis et al. 2012) to model the rel- 
ativistic effects in the spectra. Despite the existence of newer 
relativistic models that includes the black hole spin (as op- 
posed to the extent of the inner radius) as a variable parameter 
(i.e. Beckwith & Done 2004; Dovciak et al. 2004; Brenneman 
& Reynolds 2006; Dauser et al. 2010) the advantage of using 
KDBLURf on such a vast data set is it greatly expedites the 
fits and since we are not overly interested in absolute values, 
the small improvements of the newer models at the 10% level 
(Beckwith & Done 2004) does not justify their use. 

Each of the 116 remaining spectra were thus fit with a base 
model described in XSPEC as 9 

PHABS x (POWERLAW + DISKBB + KDBLURf <g> 
REFLIONX). 

To avoid unnecessary CPU cost, all spectra were fitted adopt- 
ing the same initialising values for the model parameters. A 
similar approach of using a base model with similar starting 
parameters was taken by Wilms et al. (2006) - amongst others 
- in the study of Cygnus X-l. 

Since the low-energy cutoff in the PCA of ~ 3 keV is not 
sufficient to constrain the neutral hydrogen column density, 
throughout this work we have frozen the parameter in the 
PHABS model to 7 x 10 21 atom cm~ 2 in XSPEC as the neutral 
column density is likely best modelled as being constant with 
time (Miller, Cackett & Reis 2009a) 10 . The key parameters in 
the REFLIONX model are the spectral shape of the illuminat- 
ing continuum which is set to be the same as that of the direct 
power-law and the ionization parameter of the accretion disk 
(as in Equation 1). The iron abundance was set to the solar 
value and the redshift was set to zero. 

In order to make the work more tractable, the emissivity 
profile in the blurring component was initially restricted to 
the form of a single power-law such that e(r) oc r~ q , and 
following the most recent work on the best XMM-Newton and 
BeppoSAX data of XTE J1650-500 (Walton et al. 2012) we 
have frozen the inclination of the accretion disk to 65° and 
note that this parameter is highly constrained based on optical 
light curve to be greater than 50 ± 3° and < 80° (Oroz et 
al. 2004). We also check whether the results presented below 
change if we freeze the inclination at this lower limit (50°) 
and we confirm that they do not. As is standard with such fits, 
we have frozen the outer radius to the maximum value in the 
grid of 400r g and the inner accretion disk is free to vary. 

A phenomenological combination of a LAOR emission line 

9 We have also repeated this experiment with the latest reflection model by 
Ross & Fabian (2007) (REFBHB; see also Reis et al. (2008) for a description 
of its use) which incorporates a black body component and found the results 
to be consistent in both case. We choose to carry on with the REFLIONX + 
DISKBB combination for a number of reasons; to start, this combination can 
be easily reproduced (REFBHB is not yet public) and the output parameters 
are somewhat standardised, e.g. REFLIONX has {; as a free parameter 
whereas REFBHB has the more obscure combination of hydrogen number 
density nn and the relative strength of the illuminating flux over the thermal 
emission at the surface of the disk. 

10 We used the standard BCMC cross-sections (Balucinska-Church & Mc- 
Cammon 1992) and ANGR abundances (Anders & Grevesse 1989). 
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on top of a reflection model such as PEXRIV (Magdziarz 
& Zdziarski 1995) - which does not include the iron- K a 
emission line - often provides an equally good fit to X- 
ray spectra of black holes. However, PEXRIV does not 
Compton broaden its absorption edge nor does it provide 
a physical coupling between itself and the extra LAOR 
component which can result in parameters severely lack- 
ing in physical consistency. This is exemplified in Reis 
et al. (2012) where we showed for the stellar-mass black 
hole candidate MAXI J1836-194 that the reflection com- 
ponent could be equally well modelled with a combina- 
tion of LAOR + PEXRIV; LAOR + KDBLURx PEXRIV; or 
KDBLUR x REFLIONX. The first combination unphysically re- 
quired the iron line to be coming from far within 6r g in the 
strong gravity regime, yet all other reflection features - un- 
der the formalisation of the model - appeared exempt from 
such effects. After properly correcting for relativistic effects 
in PEXRIV, the ionization parameter £ of the second model 
was nearly two orders of magnitude higher than the previ- 
ous combination as the Fe absorption edge in PEXRIV was 
no longer trying to fit the blue wing of the iron line. This 
change was also accompanied by an artificial hardening of 
the powerlaw index and a decrease in the equivalent width of 
the LAOR component from — 270 eV to — 180 eV. Correctly 
blurring the originally sharp absorption edge caused it to be- 
come smooth and more symmetric. Due to the decoupling 
between the LAOR line and PEXRIV, this smooth edge traded 
off with the LAOR line component thus decreasing the equiv- 
alent width of the latter. As such, we stress the importance of 
the imposed physical consistency in the emission line and cor- 
responding absorption edge strength afforded by REFLIONX 
and strictly enforce this in the work that follows. 

4. RESULTS AND DISCUSSION 
4.1. Evolution of the Outburst 

Figure 3 shows the evolution of all the parameters of inter- 
est in this work 11 , and Fig. 4 shows the best-fit to eight rep- 
resentative spectra roughly covering the eight periods high- 
lighted in Figs. 1,2 and described in detail in Homan et 
al. (2003). The spectra used for illustration are shown in Fig. 2 
with diamonds. The top four panels in Fig. 3 show the evo- 
lution of the extrapolated 0.1-1000kcV fluxes for the total, 
POWERLAW, REFLIONX and DISKBB components, from top 
to bottom respectively. All fluxes were obtained using CFLUX 
in XSPEC . The vertical dotted lines running through all the 
panels highlight the eight periods shown in Figs. 1,2. We see 
a clear increase in the disk flux during the first ~ 30 days 
followed by a clear flattening as it moves into the HSS. It is 
also visually apparent that the reflection flux varies relatively 
less than the power-law continuum. This will be investigated 
further in what follows. 

The following two panels show the evolution of the disk 
temperature and ionization parameter, £. Early in the outburst, 
the disk was relatively cold (< 0.5 keV) and only moderately 
ionised with £ w 200 erg cm s _1 . As the system moves into 
the HSS, the ionization increases smoothly until <~ 2 days 
before the transition to the SIS when £ sharply increases to 

1 ' We refer the reader to the work of Dunn et al. (201 1) for the evolution 
of the disk properties during the outburst. As mentioned above, despite the 
fact that the authors did not correctly account for reflection in their work, this 
plays a minor role at the energies of interest in regards to the disk properties 
and as such that work is still a valid and important reference for the evolution 
of accretion disks. 



£ w 3000 erg cm s _1 and remains at that level through the 
transition up to the end of the SIS. The disk temperature, on 
the other hand, appears to reach a relativity stable value of 
~ 0.6 — 0.65 keV early in P2, approximately half way through 
the HIS. As the system progress into the HSS, the disks be- 
comes fully ionised with £ reaching the maximum allowed 
value in the model (log £ = 4), before coming back down to 
the low hundreds towards the end of the outburst. 

The reflection fraction R shown in the third-from-bottom 
panel of Fig. 3 is here defined as a measure of the ratio be- 
tween the continuum power-law to the reflection flux emitted 
by the disk. Since a fraction of the power-law illuminating 
the accretion disk is down-scattered as it is reprocessed in the 
disk, the reflection fraction in REFLIONX is calculated by di- 
viding the extrapolated (1 eV- 1000 keV) REFLIONX flux by 
the 0. 1 - 1 000 ke V power-law flux illuminating the disk. At the 
start of the outburst, through to the end of the HIS, R increases 
smoothly between w 0.6 — 1. At the transition between the 
HIS to the SIS, R displays a sharp increase to s=s 4 where it 
remains until the beginning of the disk dominated HSS where 
the power-law effectively disappears. We also show in this 
panel the ratio of the 3-100 keV reflection to powerlaw flux. 
The behavior described above is still qualitatively the same 
and we still see a clear jump in this ratio at the transition from 
the HIS to the SIS. However, when limited to the 3-100 keV 
flux, this ratio is systematically a factor of w 1.5 less than 
the extrapolated ratio; a direct result of not accounting for the 
extra down-scattered flux at low energies. 

As a further test of both the qualitative (clear jump in R 
between HIS and SIS) as well as quantitative (change from 
R s=s 0.6 to R s=s 4 between PI and P3) behaviors found 
here for the reflection fraction, we temporally replace the RE- 
FLIONX model with a combination of LAOR plus PEXRIV and 
employ this model to the observations highlighted in Figs. 2 
and 4 for PI and P3. In using this model, we have blurred the 
PEXRIV component with the same parameters as the LAOR 
line profile Figure 5 shows the confidence range for the re- 
flection fraction (a free parameter of PEXRIV) for these two 
representative spectra. In agreement with our previous results, 
we see that in the HIS the reflection fraction is constrained to 
R = 0.58±g;"f and in the SIS it is R > 2.7 at the 90 per cent 
level of confidence (A% 2 = 2.71). 

The fact that £ is maximal in the HSS (P4) despite the fact 
that this is when the irradiation appears to be at its lowest level 
(second panel from the top) can be explained with a number of 
scenarios. We showed in § 1.1 that in stellar mass black holes 
the intrinsic hot disk can result in significant thermal ioniza- 
tion (Ross & Fabian 2007) which will be strongest in the disk 
dominated HSS. Thus, in this scenario, the high ionization 
measured is also in part due to thermal ionization. If this ther- 
mal component is the sole source of ionisation in the HSS, 
R would go to zero, consistent with the few points present in 
the HSS shown in Fig. 3. A further possibility is that the disk 
is indeed highly photo-ionised as a result of strong focusing 
of the coronal photons onto the disk. This would significantly 
remove the number of hard-photons escaping the system (thus 
explaining the second panel from the top) and cause the disk 
to be highly ionised. Indeed, observations of disk winds orig- 
inating in the HSS of various BHBs consistently show winds 
having £ <~ 10 4 erg cm s _1 (e.g. Miller et al. 2008a; Neilsen 
& Lee 2009; King et al. 2012; Ponti et al. 2012). 
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FIG. 3. — Top four panels: Extrapolated 0. 1-1000 kcV flux evolution for the total, powerlaw, reflection and disk components. The errors in the fluxes were 
assessed using CFLUX but are omitted here for clarity (but see Fig. 7 for an indication of the magnitude of the errors in each state). These errors were propagated 
in the derivation of the reflection fraction R (eighth panel from the top). We also show in the eighth panel (R) the ratio of the 3-100 keV reflection to powerlaw 
fluxes. A jump is in this ratio seen at the same position as that for the extrapolated fluxes. The bottom panel shows the reduced-^ 2 for 114 D.o.F. The vertical 
dashed lines show the transition from the rising LHS, P1-P6 and back into the falling LHS. 
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FIG. 5. — Goodness-of-fit versus reflection fraction for the two represen- 
tative spectra describing the HIS-P1 (black) and SIS-P3 (red). The spectra 
used refers to those highlighted in Figs. 2 and 4 and the REFLIONX com- 
ponent inherent in the base model has been replaced with a combination of 
PEXRIV together with LAOR. The solid blue horizontal line shows the 
90 per cent confidence range. 



4.2. Disk Emissivity 

In Fig. 6 (top-left) we show the distribution of the reduced 
X 2 (for 115 degrees of freedom) from all observations assum- 
ing the simple Newtonian 'lamp-post' like geometry in which 
the the emissivity profile follows a q = 3 power-law (light 
blue), as well as after relaxing this assumption (red). In both 
cases there is a clear peak at reduced-;^ 2 = 1, however this 
peak is much more distinct upon relaxing the Newtonian ap- 
proximation. The Newtonian approximation naturally does 
not take into account the effects of general relativity that will 
be experienced by the emission from the corona and accretion 
disk near the black hole. Relativistic effects (strong gravity as 
well as relativistic time dilation) acts to steepen the emissivity 
in the inner regions of the disk. The top-right panel of Fig. 6 
shows the distribution of the emissivity index (red) as well as 
the 90 per cent lower limit in its value (light blue). We refer 
the reader to the work of Miniutti & Fabian (2004); Wilkins 
& Fabian (2011, 2012); Fabian et al. (2012a) and references 
therein for a detailed study examining non-Newtonian values 
for the emissivity index, but note here that steep emissivities 
similar to those found here for the HIS and SIS are a natural 
and unavoidable consequence of strong gravity. 
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FIG. 6. — Top-left: Distribution of reduced x 2 f° r model with a Newtonian emissivity profile (q = 3) and a steeper values > 3. In both cases we see a peak at 
reduced \ 2 = 1> however this is much clearer after allowing the emissivity to vary beyond its Newtonian value. Top-right: Distribution of emissivity parameter 
(red) together with their 90% lower limit. Bottom: Evolution of emissivity parameter (red) together with the HEXTE-A count rate (green). The emissivity index 
was constrained to 3 < q < 10. 



The bottom panel of Fig. 6 is used here as a simple illus- 
tration of the evolution in q as well as the count-rate in both 
the PCA and HEXTE data. It is clear that q can only be con- 
strained when the PCA data is at its highest level as this con- 
straint does not come from energies > 25 keV. At the end 
of the outburst, when the PCA signal-to-noise level drops sig- 
nificantly, the data cannot differentiate between a Newtonian 
q = 3 and a steeper value. 

Following the recipe provided for Cygnus X-l by Fabian 
et al. (2012a) in dealing with sources where the spin is ex- 
pected to be high (as is the case for XTE J 1650-500), we 
have repeated our fits with a double emissivity profile such 
that within a break radius (initially frozen at 4r g but later al- 
lowed to vary) the emissivity is > 3 and beyond it is frozen 
at 3. The initial value of 4r g for the break radius was chosen 
based upon the value for Cygnus X-l (Fabian et al. 2012a). 
We find that, as long as the emissivity is not fixed at q = 3, 
the quality of the fits, and distribution of reduced % 2 are simi- 
lar to that of a single, unbroken emissivity, and we proceed by 
using this single power-law emissivity profile as our standard 
but emphasise that the results presented here do not change if 
we employ a broken emissivity profile. This is very likely to 
be due to the comparatively low spectral resolution afforded 
by RXTE which does not reflect the subtle changes in the re- 
flection profile in a similar manner as XMM-Newton or Suzaku 
observations. 

4.3. Light-Bending and General Relativity 

Hints of the expected effects of light bending, as described 
in the introduction, can be seen in the top four panels in Fig. 3. 
Most important, is the apparent constancy of the reflection 
flux in comparison with that of the direct power-law early in 
the outburst. We investigate this further in Fig. 7. The left 
panel shows the flux-flux relation through the whole outburst 
with the various spectral states shown in different colors. The 
top-right panel is a close up of the period covering the HIS 
and SIS during the first ~ 30 days of the outburst 12 . Figure 7 

12 Excluding the first 3 days when XTE J1650-500 was in the rising LHS 
(see Fig. 1). 



is remarkably similar to figure 3 of Rossi et al. (2005), where 
the authors used the flux in the iron line as a proxy for the 
total reflection in XTE J 1650-500. We have superimposed in 
this figure the expectation from the light bending model for 
a compact corona varying in height from l-20r g with a disk 
having an inclination of 60°, from Miniutti & Fabian (2004, 
see their Figure 2). The model reproduces extremely well the 
broad shape of the relation. Finally, the bottom right panel 
shows this behavior when the non-extrapolated, 3-100 keV 
fluxes are used instead. We again see that qualitatively, the 
behavior is the same as above. 

As discussed in the introduction, the light bending model 
of Miniutti et al. (2004) predicts the existence of semi-distinct 
regimes in this flux-flux relation. When the corona is located 
at a height of ~ 10r g the model predicts a flattening in the 
relation similar to that observed for the HIS (both PI and P2). 
The fluxes in this hard intermediate state are clearly corre- 
lated, and a Spearman's rank correlation test gives a coeffi- 
cient of p = —0.799 corresponding to a 1.7 x 10~ 6 chance 
of a false correlation (Table 1). The slope of this relation is 
—0.24 ± 0.02 (standard error, s.d.) and this linear fit is shown 
in the top-right panel of Fig. 7 as a solid black line. 

As the location of the corona reaches the more extreme en- 
vironment within a few r g from the black hole, the model 
predicts a steep, positive linear relation between the reflection 
and power-law flux similar to that seen in the SIS, although 
there is large scatter. A Spearman's rank correlation test here 
gives p = 0.719, with a false correlation probability of just 
2.4 x 10~ 6 . The slope of this relationship is greater than unity, 
with a best-fit value of 2.7±0. 1 (shown as a further solid black 
line in Fig. 7) . Note that this is highly inconsistent with the 
expectation for a static Newtonian corona with intrinsically 
varying luminosity, where the slope should be unity across the 
entire flux range. It is the combination of a slope ^> 1 at low 
powerlaw flux together with a near-flattening at higher fluxes 
that provides evidence for relativistic effects in this case. 

It has been suggested (e.g. Ballantyne et al. 2003, 2011) 
that changes in the ionization of the inner regions of the accre- 
tion disk can give rise to changes the reflection spectrum that 
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FIG. 7. — Left: Flux-Flux relation during the outburst. The dashed line shows the one-to-one relation. Top-right: Close up of the HIS (blue and cyan) and SIS 
(green). The solid black lines shows the best linear fit for each state having a slope of 2.7 ±0.1 and —0.24 ± 0.02 in the SIS and HIS respectively. The solid 
red curve shows the expected flux-flux relation under the light bending model of Miniutti & Fabian (2004, see their Figure 2) for a system with 60° inclination, 
somewhat similar to XTE J1650-500, with a corona varying in height from l-20r g . Note that this is not a fit and has been rescaled from the original. The solid 
red dots indicate a corona height of 1 r g and 5r g from left to right. Bottom-right: Similar but for fluxes between 3-100 keV. 



TABLE 1 

Summary of Spearman's Rank Correlations and Partial Correlation tests on a combination of various model parameters. 





HIS (PI + P2) 


SIS (P3) 






Spearman's 


Partial 


Spearman's 




Partial 






rank-order 


correlation 


rank-order 


correlation 


Parameter 1 


Parameter 2 


P 


p-value 


P 


p-value 


P 


p-value 


P 


p-value 




-^powcrlaw 


-0.953 


2.0 X 10" 13 


-0.473 


0.014 


0.319 


0.071 


0.092 


0.621 


S 


^disk 


0.956 


9.3 x 10~ 14 


0.515 


0.006 


0.006 


0.972 


0.082 


0.656 


-fdisk 


-^powcrlaw 


-0.955 


1.1 x 10~ 13 


-0.427 


0.0031 


-0.310 


0.079 


-0.334 


0.056 




-FrcflionX 


0.742 


2.2 x 10- s 


-0.140 


0.516 


0.382 


0.028 


0.218 


0.230 


^disk 


-FrcflionX 


0.766 


7.8 x 10~ 6 


0.088 


0.685 


-0.116 


0.520 


0.140 


0.446 


-^powcrlaw 


-FrcflionX 


-0.799 


1.7 x 10" 6 


-0.373 


0.066 


0.719 


2.4 x 10" 6 


0.683 


4.6 x 10~ 7 



Notes:- Spearman's rank correlations and Partial Correlation test were made for combinations of the reflection, power-law and disk fluxes as well as the 
ionization parameter. The Partial Correlation test measures the degree of associations between the two parameters listed on the first two columns whilst 
controlling for the remaining two parameters. The Spearman's coefficient p is a measure of the degree of correlation with +1 or -1 indicating a perfect monotone 

function and a lack of correlation. 



can mimic somewhat this flat behavior. In order to robustly 
assess the strength of the correlations seen in Fig. 7 in both 
the HIS and SIS, we have performed a number of correlation 
tests which are summarised in Table. 1. We performed, for 
both HIS and SIS, Spearman rank-order tests for a combina- 
tion of four parameters (power-law, reflection and disk fluxes 
as well as the disk ionization) as well as Partial Correlation 
tests (PCT) for two parameters while controlling for the third 
and fourth variable. The PCT is of particular importance for 
our purpose as it removes any potential association of the ion- 
ization parameter (or any other potential source of unwanted 
correlation) in the flux-flux relations shown in Fig. 7. 

From the Spearman rank-order tests performed in the HIS 
(Table 1), it would initially appear that all four variables are 



strongly correlated with one another in some way, as all com- 
binations display \p\ > 0.7. However, after performing the 
partial correlation test for all combinations we see that for 
two of the previous strong correlations (£ — F Ic fn on x and 
^disk — ^reflionx) were in fact driven by the mutual depen- 
dence of these parameters on F powc \ aw . These tests clearly in- 
dicate that the reflection flux in both states is better correlated 
with the power-law flux than with ionization and the slope of 
the correlations (mildly negative in the HIS and strongly pos- 
itive in the SIS) are highly indicative of gravitational light- 
bending in the General Relativistic regime. 

The behavior seen here is consistent with a drop in the 
height of the corona during the hard-intermediate phase (PI 
and P2) followed by intrinsic variations in its luminosity by 
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a factor of a few during the soft-intermediate state (P3). Fol- 
lowing the disk-dominated soft-state (P4), the height of the 
corona increases again (P5 and P6) and the outburst finishes 
with the intrinsic power dropping as the source goes back into 
the I.I IS. 

4.4. R — T relation 

A strong correlation has been shown to exist between the 
amplitude of the reflected component {R) and the photon in- 
dex (T) of the Comptonized spectrum in XRBs in the hard 
state (Ueda et al. 1994). This R — T relation has since been 
robustly tested by a number of authors (e.g. Gilfanov et al. 
1999; Zdziarski et al. 1999, 2003; Nowak et al. 2002; Ibragi- 
mov et al. 2005) and it is now thought to also apply to Seyferts 
and radio galaxies (e.g. Zdziarski et al. 1999), further cement- 
ing the similarities in the coronal properties at all scales. If 
this relation indeed turns out to be real (and evidence attests 
to this; but see Molina et al. 2009) then this is could be telling 
us about the feedback between the hot corona and cold gas in 
an accretion disk. 

We briefly investigate this relationship for XTE J 1650-500 
in Fig. 8 where we show using a combination of a 
phenomenological model for the Comptonisation with a 
physically-consistent prescription for the reflection compo- 
nent, that the R — V correlation is real but only persists during 
the low-hard and early hard-intermediate states. 

The presence of this relation early in the outburst and af- 
ter the system returns to the LHS strongly points towards a 
feedback process between the soft photons in the disk and the 
corona. There are a number of theoretical interpretations for 
the presence of this correlation (Poutanen et al. 1997; Gil- 
fanov et al. 1999, 2000; Beloborodov 1999b; Malzac et al. 
2001) with the two leading contenders often described as the 
disk-truncation and dynamic corona model Done see 2002; 
Beloborodov see 1999a, for a detailed study of these models. 
To summarise, in the former, an increase in the reflection frac- 
tion is caused by the accretion disk penetrating deeper into a 
central hot corona thus receiving more illuminating hard pho- 
tons. The presence of the disk in return offers more soft pho- 
tons consequently cooling the plasma. For a purely thermal 
distribution of electrons, the greater the number of soft seed 
photons, the softer the power-law spectra. The latter model 
by Beloborodov (1999b) invokes bulk-motion of the corona 
above the accretion disk. If the corona is outflowing with 
mildly relativistic velocities, this would reduce the amount 
of hard photons hitting the disk, which in turn reduces R and 
subsequently softens T as few reprocessed soft photons reach 
the outflow. Recent evidence for the coronal plasma ejection 
model has come from a strong positive correlation between 
reflected X-ray flux and radio flux in the black hole binary 
Cygnus X-l (Miller et al. 2012). 

We will show in § 4.6 that all evidence points toward the 
disk radius remaining stable during the HIS to SIS transi- 
tion (see Figs. 1 1,12). This constancy effectively rules out the 
"disk-truncation" explanation for the R — T correlation. In- 
stead, one may hypothesise whether the "outflowing corona" 
and light bending can be combined to explain the behavior so 
far detailed for XTE J 1650-500. In the previous section we 
showed that the flux-flux behavior (Fig. 7) could be explained 
if early in the outburst (during the HIS) the corona was lo- 
cated relatively far (~ 10r g ) from the black hole and thus 
behaved according to Regime 2 of Miniutti & Fabian (2004, 
see also §1.2). As the system evolved, the corona collapsed to 
a few r g and began to experience a higher level of light bend- 
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FIG. 8. — The well know R-gamma relation is clearly present in the LHS 
(with a few outliers) and during most of the HIS. However, just prior to the 
transition to the SIS and thereafter, this relation does not hold. We discuss a 
potential explanation for this behavior in § 4.4. For clarity, we have removed 
the errorbars for the points in the LHS. 

ing towards the disk (regime 1). In this scenario, a potential 
gradient in the outflow velocity of the corona as a function 
of height could also explain the behavior seen in Fig. 8. i.e. 
as the corona collapses from a large height (large outflow ve- 
locity; low R and hard T) the outflow velocity decreases (R 
increases and T softens) until it becomes effectively static and 
the system transitions into the SIS. We will expand upon this 
possible scenario in what follows and summarise our ideas in 
§5. 

4.5. QPOs and Spectral States: A collapsing corona 

Throughout this work, we followed the selection made by 
Homan et al. (2003) which roughly separates the outburst 
into six periods coinciding with significant changes in both 
the hardness intensity diagram (Fig. 2) as well as in the 
shape of their power spectra 13 As highlighted in the introduc- 
tion (§1.3), those authors demonstrate the presence of high- 
frequency (HF) variability in XTE J 1650-500, together with 
a HFQPO which was shown to evolve in both frequency and 
coherence during the outburst. The highest frequency which 
was reliably measured was at ~ 250Hz in the SIS, with the 
frequency being much lower (~ 50Hz) at the onset of the out- 
burst. 

In Fig. 9 (top), we show the presence of a strong (Pearson's 
r = 0.938) positive relation between the reflection fraction 
and quality factor of the QPO (Q — R relation). In order to 
create this figure, we have averaged the values of R shown in 
Fig. 3 for each of the periods in question and added a sym- 
metric error of ±1 in the quality factor during Period 5 as this 
error was not provided by Homan et al. (2003). The bottom 
panels shows the frequency of the HF QPO as a function of R 
(/ — R relation). 

QPOs are notoriously difficult to explain and it is not our 
purpose to provide a quantitative description of this phe- 
nomenon. However, it is worth stressing that most models 
(e.g. Nowak et al. 1997; Cui et al. 1998; Psaltis et al. 1999; 
Stella et al. 1999; McKinney et al. 2012) strongly link the 

13 In keeping with that work, the HIS is divided into two periods (PI and 
P2). Here, we also add two extra periods which we have denoted as the LHS- 
rising and falling. 
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FIG. 9. — Top: HF QPO coherence as a function of reflection fraction. The 
relation is well described by a linear function similar to that shown in the fig- 
ure. Bottom: QPO frequency as a function of reflection fraction. The dashed 
line shows a relation, /(Hz) = (102 ± 2) X ln[i? X (3.07 ± 0.05)], that fits 
these data. These figures show a clear link between the coherence/frequency 
of the QPO and the reflection fraction which we subsequently interpret as 
being linked to the size/position of the corona (see § 4.5). 

origin of QPOs with orbits and/or resonances in the inner ac- 
cretion disk close to the black hole. Current models cannot 
fully explain, in a physical manner, the range in coherence 
observed in various systems nor the manner in which the fre- 
quencies change with states. 

To explain the range in coherence observed in accreting 
neutron stars, Barret et al. (2006, 2007) devised a toy model 
in which the changes in Q are driven by changes in the scale 
height of the disk. A small scale height gives rise to high co- 
herence and vice versa. Expanding on this idea, it appears 
that at least for XTE J 1650-500, it is physical changes in the 
radius/size of the corona that give rise to changes in both qual- 
ity factor Q and QPO frequency. To illustrate this hypothesis, 
consider Fig. 9 (bottom) together with our interpretation of the 
behavior displayed in Fig. 7 (§ 4.3). Early in the outburst the 
frequency of the HF QPO appeared at ~ 55Hz. The Keplerian 
frequency at a given radius is /(Hz) sa 3.2 x 10 4 AI~ 1 r~ 3 / 2 , 
where M and r are in units of Solar mass and Gravitational 
radius respectively. Hence, during the brief PI period, the HF 
QPO frequency is close to the orbital frequency at ~ 28r g as- 
suming a 4 M black hole and potentially moves to ~ 15r g in 
the second half of the HIS. As the outburst continues and the 
corona continues to collapse, it is plausible that the frequency 
continues to increase (corresponding to ~ 10r g in the SIS) 
eventually approaching a value that should be consistent with 
the Keplerian frequency at the ISCO. The continued decrease 
in the size of the corona gives rise to the increase in coher- 
ence. In this scenario, the frequency of the QPO should relate 




2 3 
Reflection fraction (Ji) 

FIG. 10. — Radio flux density as a function of reflection fraction. The 
arrows show the direction of the outburst. There appears to be two distinct 
branches where the reflection fraction is either constant at R ~ 0.5 or at R ~ 
4. These two branches corresponds to the LHS/HIS and SIS respectively. 
Radio data obtained from Corbel et al. (2004). Corbel et al. (2004). 

to the size of the corona and thus would naturally increase as 
the corona collapses. 

The relation between the coherence of the HF QPOs and 
the reflection fraction also leads to the interesting prediction 
that HF QPOs should only be observed when R > 0.4 when 
the coherence is Q > 0. This results is consistent with obser- 
vations, with no HF QPO ever having been found in the LHS 
where R < 1. 

4.6. Radio (jet) Emission and Reflection Fraction 

Corbel et al. (2004) presented a comprehensive analy- 
ses of the radio emission observed during the outburst of 
XTE J 1650-500. In that work, the authors suggest that the 
transition between the HIS and SIS 14 is associated with a mas- 
sive radio ejection event. The observations in the LHS was 
found to be consistent with the presence of a steady, compact 
jet, as is often seen in the LHS of black hole binaries (see 
e.g. Fender 2001; Fender et al. 2009, and references therein). 
This potential ejection event during the HIS-SIS transition 
coincides with the time where we see a sharp jump in the 
reflection-fraction. In Fig. 10 we show the radio flux density 15 
vs the reflection fraction calculated herein. At the time of the 
steady compact jet (in the LHS & HIS), the radio flux density 
increased by a factor of ~ 5 with no statistically significant 
change in R. However, immediately following the radio ejec- 
tion in the SIS, the reflection fraction increases dramatically 
resulting in a bi-modality in the radio-flux density-reflection 
fraction plane. This suggests an intimate link between the jet 
ejection site and the collapsing corona. Indeed, Beloborodov 
(1999b) predicts a link between radio jets and reflected flux, 
which was also previously seen in Cygnus X-l (Miller et al. 
2012). 

At later stages in the outburst (HIS-P6, falling LHS), the 
measured radio flux density is likely a dominated by emission 
from the zone where the ejected plasma interacts with the am- 
bient ISM surrounding the binary system. Unfortunately, due 

14 Referred to as the intermediate and steep powerlaw states respectively 
in Corbel et al. (2004). 

15 Obtained directly from Table 1 of Corbel et al. (2004). We are using 
the flux densities at 4800 MHz for all observations, except for the rising-LHS 
where this was not available. In this case, we proceeded by averaging the 
values presented for 1384 and 2496 MHz. 
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FIG. 1 1. — Left: The inner radius during the outburst are shown in both red and blue, where the blue points indicate where the uncertainty in the radius 
measurement is < 50% of its value. The dashed lines show the evolution of the PCA (grey) and HEXTE-A (green) count rate in a similar manner to Fig. 6 
(bottom). Right: Distribution of inner radii having errors < 50% its value. In red we show the standard normal distribution with a mean radius of 1.65r g and 
standard deviation of 0.08rv. 




Link Number Link Number '« Number 

FIG. 12. — MCMC results for the inner radius (a proxy for spin) obtained from the simultaneous fit to the first ten spectra in the HIS. Left: Figure tracing 50 
out of a total of 170 "walkers" during their random walk. The figure shows that the various chains converge pretty quickly indicating efficient sampling. The 
inset shows a close up of the first 20,000. steps. Right: Full MCMC containing all 170 walkers, after having the first 5,000 elements "burnt-in". For clarity, we 
only show every 1000th element of the chain. 



to the low spatial resolution of the radio observations, this 
emission is not resolved from that due to any reformation of 
the steady jet close to the black hole. 

4.7. Disk ( inner) Radius and State Transition 

As has been discussed throughout this paper, a popular ex- 
planation for state transitions is a radial variation in the ex- 
tent of the accretion disk. This model has been highly suc- 
cessful in part due to its flexibility and ease in which it can 
explain the "weak" reflection fraction (R < 1) often found 
in the LHS. However, over the past few years it has consis- 
tently been shown that in the luminous phases of the LHS - 
at least above ~ 1 x 10~ 3 L^dd 16 - the disk does not ap- 
pear to be truncated away from the radius of the ISCO (Miller 
et al. 2006b,a, 2008b; Reis et al. 2009c, 2010; Reynolds & 
Miller 2011; Walton et al. 2012). Figure 11 (left) shows that 
the present work cannot statistically rule out a disk being trun- 
cated late in the outburst, once the system has returned to the 
LHS, however it strongly constrains this radius to ~ 1.65r g 

16 Contrast this with the broadband analyses of GX 339-4 presented by 
Tomsick et al. (2008) where the authors find clear evidence for the reces- 
sion of the accretion disk beginning only at Eddington luminosities below 
~ xlO" 3 L E dd- 



during the transition from the HIS to the SIS (right), thus rul- 
ing out radial variation as the cause for the HIS-SIS transition. 

Where we have not been able to constrain the radius, this 
has largely been due to the data quality (the LHS is inherently 
less luminous) as well as the fact that reflection is intrinsi- 
cally weaker in the LHS. The strongest reflection features are 
expected in the intermediate states where the disk receives 
a larger fraction of the hard X-ray emission (see e.g. Hiem- 
stra et al. 201 1). Note also that despite the comparatively low 
spectral resolution afforded by the RXTE-PCA (18% FWHM 
energy resolution at 5.9 keV) Wilms et al. (2006) showed 
that this instrument can indeed resolve line-widths down to 
a least a ~ 0.3 keV. Higher resolution observations of 
XTE J 1650-500 early in the outburst showed that, when mod- 
elled with a GAUSSIAN, the FeK a emission line is consistent 
with having a width of a ~ 1.1 keV (Table 4 in Walton et al. 
2012) 17 . 

Assuming the stable radius shown in Fig. 1 1 for the two in- 
termediate states is indeed the radius of the ISCO, we find, 
using the relationship between ISCO and black hole spin 

17 The original analyses of this XMM-Newton dataset performed by Miller 
et al. (2002) included an extra smeared edge component at ~ 6.8 keV which 
resulted in the GAUSSIAN having a width of only ~ 250 eV. 
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of Bardeen et al. (1972), a dimensionless spin parameter of 
0.977^q'qqy consistent with the value found in detailed anal- 
yses of single, high quality data obtained with XMM-Newton 
(0.84 < a < 0.98; Walton et al. 2012) or BeppoSAX (a > 
0.93; Miniutti et al. 2004 18 ). 

As a test of the robustness of this result, we have performed 
a joint fit to the first ten observations in the HIS. We used the 
same base model as before with each individual observation 
having their own set of parameters - disk emissivity index, 
temperature, normalisation and parameter, power-law index, 
as well as the normalisation of the powerlaw and reflection 
component. However, this time we forced the inner radius in 
the various observations to be a global parameter thus assum- 
ing a constant value. 

This simultaneous fit contains a total of 81 free param- 
eters 19 and with this comes a high chance of mistaking 
a local minima in % 2 space for the global best fit. In 
order to address these limitations, we proceeded by min- 
imised the fit using standard x 2 fitting techniques within 
XSPEC until a reasonable fit was produced (x 2 /v < 2) 
at which point we halted the minimisation 20 and proceeded 
with Monte Carlo Markov Chain (MCMC) analysis. We 
employed the MCMC procedure described in Foreman- 
Mackey et al. (2012) (code found at http://danfm.ca/emcee/) 
and implemented in the XSPEC spectral fitting package 
by Jeremy Sanders (XSPEC implementation described in 
https://github.eom/j eremys anders/xspec_emcee) . MCMC 
techniques have been successfully used to address similar 
problems in constraining the black hole spin of NGC 3783 
(Reynolds et al. 2012) as well as in modelling the kinematics 
of the microquasor XTE J 1550-564 (Steiner & McClintock 
2012). 

We added a 5% random perturbation to all the parameters 
in the fit described above, and increased the value of the in- 
ner radius from the starting value of <~ l-6r g to 2.5r g in or- 
der to guarantee that the chain could freely converge to the 
global minimum. We used a total of 170 "walkers", each it- 
erated ("walking") 10,000 times. Figure 12 (left) shows the 
evolution of the walk in inner radius for 50 randomly selected 
walkers. It is clear that the walkers converge to the same value 
efficiently. Nonetheless, in order to be conservative we have 
ignored ("burned-in") the first 5000 elements of each chain 
and show on the right the full MCMC chain for the inner ra- 
dius which is clearly well behaved with a peak distribution at 
1.66 ± 0.005r g (s.d.), in excellent agreement with the results 
in Fig. 11. 

5. SUMMARY 

In this work, we try to take a broad systematic approach to 
not only the data reduction but also in the manner in which 
the spectra are fit during all spectral states observed during 
the outburst. We have presented a number of empirical results 
inferred from our a priori assumption that the observed spec- 
tra and their evolution is a consequence of variation in three 
separates emission components - the power-law continuum, 
thermal disk and reprocessed reflection. Although we have 
tried to convey a fully consistent, albeit qualitative picture 

18 Spin converted from the lower limit on the inner radius of ~ 2.1r g . 

19 The sheer number of free parameters and computational time required 
to do x 2 fitting as well as the MCMC analyses described in what follows 
drove the need to constrain this analyses to only 10 observations as opposed 
to all 116. 

20 The actual quality of the fit at this time was \ 2 / u = 2059.7/1149. 



of the spectral evolution of XTE J 1650-500, it is very diffi- 
cult to obtain unique interpretation of the results and nearly 
impossible to generalise this to all systems. In a forthcom- 
ing publication, we will apply similar techniques to a larger 
set of objects at which point we hope to be able to make a 
stronger statement regarding the global population of stellar 
mass black hole binaries and possibly AGN. 

For now, we summarise the main results of the presented 
spectral analyses as follows: 

1 . The outburst is well characterised by a model consisting 
of a hard X-ray continuum, a thermal disk component 
and reprocessed emission (reflection). 

2. The emissivity profile of the disk is not well char- 
acterised by a simple Newtonian approximation and, 
where this value can be constrained early in the out- 
burst, it is steeper than the Newtonian value q = 3. 

3. The Flection - -Fpowcriaw plane for the periods cover- 
ing the hard-intermediate to the soft-intermediate spec- 
tral states displays two distinct behaviors. Early in the 
outburst, during the HIS, this plane is nearly flat (slope 
of —0.24 ± 0.02) with the powerlaw flux varying by 
<~ 5x and the corresponding reflection flux by <~ 1.5 x 
During the SIS that followed, this behavior become dis- 
tinctively different with a F ro fl cction - F powcrlaw slope 
of 2.7 ±0.1. 

4. The nearly-flat behavior of the inflection - -Fpowcriaw 
plane seen for the HIS cannot be explained away by 
variations in the ionization of the accretion disk and we 
propose that the most likely explanation is that of the 
the light bending model of Miniutti & Fabian (2004). 

5. The HIS-SIS transition is accompanied by a sharp in- 
crease in the reflection fraction, which we interpret as a 
sudden collapse of the corona as the system approaches 
the thermal state. The collapsed corona now experi- 
ences stronger effects of light-bending due to its prox- 
imity to the black hole. The radiation from the corona is 
focused towards the disk thus systematically increasing 
the reflection and decreasing the continuum. 

6. We confirm the R — T correlation during the low-hard 
and early hard-intermediate states but find that this re- 
lation does not hold once the system has transited into 
the SIS. 

7. We find a strong linear correlation between the reflec- 
tion fraction and the coherence of the high frequency 
QPOs. We also find a correlation between the fre- 
quency of the QPO and R which is well explained with 
a simple log-linear relation. 

8. We have presented a scenario within the collapsing 
corona toy model where the increase in the coherence 
of QPO with increasing reflection fraction is a conse- 
quence of decreasing emitting region as the corona col- 
lapse to regions closer to the black hole where light 
bending significantly increases R. Similarly, as the 
corona collapses, the increase in the QPO frequency 
could be associated with a decrease in the radial extent 
of the corona. 
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9. In the low-hard and hard-intermediate states, the radio 
flux density varies by a factor of <~ 5 with no change in 
the reflection fraction. The "cannonball" radio ejection 
associated with the HIS-SIS transition is accompanied 
by a sharp increase in the reflection fraction. 

10. We show that the HIS-SIS transition is not due to vari- 
ations in the inner radius of the accretion disk which is 
found to be stable at 1.65 ± 0.08r g . Assuming this is 
the radius of the innermost stable circular orbit, we find 
a spin of a > 0.96 in excellent agreement with values 
measured in other works. 
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